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We provide a theoretical analysis for the intermediate scattering function typically measured in depolarized 
dynamic light scattering experiments. We calculate the field autocorrelation function t) in dependence 

on the wave vector Q and the time t explicitly in a vertical-horizontal scattering geometry for differently 
shaped solids of revolution. The shape of prolate cylinders, spherocylinders, spindles, and double cones with 
variable aspect ratio is expanded in rotational invariants fim{'>')- By Fourier transform of these expansion 
coefficients, a formal multipole expansion of the scattering function is obtained, which is used to calculate 
the weighting coefficients appearing in the depolarized scattering function. In addition to translational and 
rotational diffusion, especially the translational-rotational coupling of shape-anisotropic objects is considered. 
From the short-time behavior of the intermediate scattering function, the first cumulants UiQ) are calculated. 
In a depolarized scattering experiment, they deviate from the simple proportionality to Q^. The coefficients 
fimiQ) strongly depend on the geometry and aspect ratio of the particles. The time dependence, in addition, 
is governed by the translational and rotational diffusion tensors, which are calculated by means of bead models 
for differently shaped particles in dependence on their aspect ratio. Therefore, our analysis shows how details 
of the particle shape—beyond their aspect ratio—can be determined by a precise scattering experiment. 
This is of high relevance in understanding smart materials which involve suspensions of anisotropic colloidal 
particles. 
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I. INTRODUCTION 


Although the thermal motion of small suspended particles is connected to the name of Robert Brown, who observed 
a random motion of pollen in water with an optical microscope in 1827, this effect was for the first time reported by 
Jan Ingenhousz in 1785. More than 40 years earlier, Ingenhousz described the stochastic motion of inorganic charcoal 
dust on the surface of alcohol, already disproving Brown’s hypothesis that this motion is caused by the vitality of 
pollen. The physical reason for this thermal motion was elucidated by Einstein in the annus mirabilis of physics, 
1905.13 Already one year later, he addressed the rotational diffusion of spherical particles.13 

While the translational diffusion of colloidal particles has been thoroughly explored by experiments and theory,l331 
rotational diffusion is much less investigated. The early starting point for the lasting interest in the diffusive mo¬ 
tion of anisotropic particles is the well-known work of Perrin.13 Experimental access to rotational diffusion is possible 
via dynamic light scattering employing laser lightPti^ or coherent X-ray^i3 as a probe. With near-held depolarized 
dynamic light scattering, in a heterodyne experiment, very small scattering vectors can be accessed.^ Eurthermore, 
birefringenc^i3 and NMR experiments^ as well as computer simulationJi^lini are used to investigate rotational diffu¬ 
sion. For particles with a size comparable to optical wavelengths, the rotational and translational motion can directly 
be observed by means of confocal video microscopy.l^ As opposed to scattering methods, which immediately provide 
information on an ensemble average, video microscopy requires tracking a large number of trajectories for a statistical 
average. Furthermore, these methods are limited to comparatively large and slowly diffusing objects. 

While a variety of well-dehned spherical particles have been studied, only few model systems of well-dehned 
anisotropic particles are available. Systems with tunable aspect ratio are even rarer. As model systems for pro¬ 
late particles, e.g., goethite,^^!}!!!] ]^gj^g^|.j.(.g| 26 H 28 ] pgjyj^gj. particles,^ and viruses like tobacco mosaic viruse^^^ESl or 
fd viru ses^^ BSl j^g^.^g been described. With regard to oblate particles, gibbsite platelets are available as a model 

system .1^3111 

Due to the orientation as an additional degree of freedom, however, the complexity of the phase diagram of such 
anisotropic particles is enhanced as compared to spherical particles since at least the interaction by the excluded 
volume is dir ection-dependent. Both prolate and oblate particles with large anisotropy are mesoscale analogs of liquid 
crystal^^^ljoi therefore exhibit the same variety in phase behavior. 

Anisotropic colloidal particles are important for many complex fluids in technical applications, e.g., suspended 
pigments in paints, clay particles in slurries used for the preparation of ceramics, and even in food engineering.SH 
Ferrofluid^^ and ferrogel^^^EH (for a recent review see Ref. with magnetic particles provide further prominent 
examples of smart materials where anisotropic particles are involved. Both translational and rotational friction of 
those particles influence the macroscopic rheological behavior important for the processing of such fluids. Since the 
size of colloidal particles is similar to that of biological structures such as viruses or bacteria, their microrheology is 
also relevant to biological processes. Shape-anisotropic particles, even when consisting of amorphous materials with 
isotropic physical properties, exhibit direction-dependent dielectric and magnetic properties and herewith related 
optical properties. Hence, the rotational motion of such particles affects their scattered intensity. 

The translational self-diffusion can be addressed experimentally by quasielastic scattering experiments using longitu¬ 
dinal and transversal waves without polarization analysis. For the experimental observation of rotational self-diffusion 
by means of scattering methods, however, transversal waves with polarization analysis are mandatory as a probe. The 
rotational motion can be elegantly probed using an incident wave with a polarization perpendicular to that of the 
scattered wave. If the polarization of the incident beam is perpendicular to the scattering vector, the geometry of 
such a setup is called vertical-horizontal. 

Let us consider an ensemble of freely rotating particles completely aligned along an initial unit vector Dg at time 
t = 0. This ensemble is described by a sharp orientational distribution function (odf) p(r2) = 5{d — Dg) with fl 
denoting a unit vector on the sphere. Clearly, this orientational distribution function becomes random with p(r2) = 
l/(47r) after a very long time. Hence, the mean angle between the initial and final orientation, 
equals tt/2, which corresponds to the angle enclosed by the crossed polarizers in a depolarized scattering experiment. 

The key dynamical quantity measured in a depolarized light scattering experiment in vertical-horizontal (VH) 
geometry is the so-called intermediate scattering function, which is related to the field autocorrelation function and 
provides access to the orientational dynamics. Even for a single Brownian particle, i.e., for a suspension at high 
dilution, the time dependence of this correlation function is nontrivial due to the rotation-translation coupling in the 
laboratory frame. A first theoretical approach to the translational-rotational coupling of cylinder-shaped objects is 
described by Aragon and PecorgP^using spheroidal basis functions. In this paper, we start from a systematic expansion 
of the particle form factors P{Q, O) into spherical harmonics. The intermediate scattering function in VH-geometry 
is then calculated by canonically averaging over rotational and translational transition probabilities obtained from 
the Smoluchowski equation for the Brownian dynamics of anisotropic particles. We give explicit expressions for the 
expansion coefficients and evaluate them for different particle shapes such as cylinders, spherocylinders, ellipsoids, 
spindles, and double cones. We thereby show that details of the particle shape—beyond the aspect raticPSl —influence 
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the static and dynamic scattering functions. Therefore, a precise scattering experiment, if evaluated within our 
framework, gives access to the details of the particle shape. This is of high importance with regard to the interpretation 
of depolarized quasielastic scattering experiments employing coherent X-rays. When colloidal suspensions are studied, 
the size of the particles is significantly larger than X-ray wavelengths. Meanwhile, highly brilliant linearly polarized 
coherent X-rays are accessible from third generation synchrotron source^^and new upcoming free electron lasers .1^ 
This paper is dedicated to the depolarized scattering of solids of revolution. Our approach can be applied to prolate 
as well as to oblate shapes. Here, at first prolate objects are considered. 


The paper is organized as follows: in Sec. 
terms of rotational invariants. In section III 


^ we calculate the field autocorrelation function in VH-geometry in 
the scattering length density of homogeneous solids of revolution is 


expanded in spherical harmonics. The static scattering functions are the Fourier transforms of these formal multipole 
expansions as described in Sec. |IV[ The method which is used to calculate the translational and rotational diffusion 
coefficients for the various considered particle shapes is presented in Sec. |V] These diffusion coefficients are needed 
to describe the time dependence of the depolarized scattering function. The resulting first cumulants are analyzed in 
Sec. |VI[ before we finally conclude in Sec. ivnl 


II. DEPOLARIZED FIELD AUTOCORRELATION FUNCTION FOR SOLIDS OF 

REVOLUTION 


Let p(r, r2,t) denote the probability density to find a particle at the time t located with its center of mass at the 
position r and with orientation D = (sin z? cos v?, sin sin i^, cos i?). The orientation O is characterized by the angles 
{z?, (/?}, where z? is the polar angle and ip the azimuthal one. 

The translational and rotational diffusive motion of an ensemble of non-interacting particles can be described via 
the Smoluchowski equatioiJ^ 



. V r -|- D^ot 



p(r,D,t) 


m 


p(r,D,t), 


( 1 ) 


where Dy and D± are the translational diffusion coefficients parallel and perpendicular to the rotation axis of the 
body of revolution, while is the rotational diffusion coefficient for Brownian reorientations of the particle. Vr is 
the spatial gradient, and the operator 


is the angular part of the Laplacian acting on the orientation D of the particle. 
In Fourier space, the Smoluchowski equation reads 


J_ ^ 

sinz? Uz? 5z? 


1 

sin^ z? 


( 2 ) 




+ 2 AUtr 



.Q + DroJ"^ p(Q,D,t) 




p{Q,n,t ), 


(3) 


where Q is the scattering vector. The translational diffusion tensor is decomposed in irreducible spherical tensors 
T(°^ and T*^^^ using the abbreviations (Dtr) = {2D± -|- D||)/3 and ADtr = Dy — D±. T*^°^ is the unit tensor in the 
first term of Eq. ([^ and T^^^ is the traceless symmetric tensor in the second term, corresponding to formal monopole 
and quadrupole contributions to the translational diffusion. The contraction of the irreducible spherical tensors with 
the unit vectors Q leads to 

(^-Q^(Ar) - I g" AAr P 2 (cos d) + piQ , t) = ^^p{Q , D, t) , (4) 

where the polar angle d corresponds to the angle enclosed between the particle axis assumed to be directed in the 
positive ^-direction and the scattering vector Q = Qfl. Since scattering methods probe the projection of translational 
diffusive displacements on the scattering vector, the angle z? is also the angle enclosed between the particle’s orientation 
and the direction of its diffusion path. A (a:) = (3a;^ — l)/2 is the second Legendre polynomial with x = cosz?. By 
means of the decomposition of the diffusion tensor in irreducible spherical tensors, the Smoluchowski equation in 
Fourier space for the translational and rotational diffusion of axially symmetric solids of revolution can formally be 
written as 

Cp{Q, n, t) = + Aot) p(g, t) = ^PiQ, t) 
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(5) 
















with = —Q'^iDtr) accounting for the isotropic or monopole part of the translational diffusion tensor, CADtr = 

— (2/3)(5^AZ)tr^2(cosi?) for its quadrupole part, and, finally, Crot = D^ot^ for the rotational diffusion. 

The solution of this parabolic partial differential equation is formally given by 

p{Q, n,t) = e^^piQ, n, = e^<°->‘e^^-t^‘eA°t*p(Q, n, , (6) 

since the operators CADt^, and C^ot commute. 

The depolarized scattering function gY^{Q,t) is determined by the product of the solutions of the three parts of 
the Smoluchowski equation. As the different contributions are statistically independent, the scattering function can 
be written as 


9riQ,t) = 9lfD,^){Q,t)gl'inJQ,t)9rUQ,t) ■ 


(7) 


Here, ^ {Q, t) accounts for the isotropic part of the translational diffusion, while g^^o^ (Q, t) gives the direction- 
dependent contribution. Finally, the effect of rotational diffusion is represented by the term la the 

following, these parts are considered separately. Since for the calculation of the quadrupole contribution to the 
translational diffusion the results of both the rotational and the isotropic translational part are needed, the following 
two subsections briefly review corresponding textbook knowledge.^ 

The depolarized field autocorrelation function is for an ergodic system given bjEi 


gX^iQ^t) 


(E*(g, ^ 0 ,0)G(O, t\no, 0 )E(Q, n, 
(E*(Q,no, 0 )E*(Q,f 2 o, 0 ))j,^ 


( 8 ) 


where the notation (.. indicates canonical averaging over the initial orientation fig the final orientation H 

of the particles. The scattered electromagnetic field expanded in spherical harmonics reads 


E(g,n(t)) = EE 

I m 


(9) 


and E*(Q, r2(t)) is its complex conjugate. The kernel G(r2,t|Oo,0) is the conditional probability density that an 
object’s orientation is fl at time t given that it was fJo at the time to = 0. This conditional probability density is the 
product of the transition probability p(Jl, 0) for a rotation from Og to fl and the probability density to find the 
particle in a statistical ensemble with the initial orientation rio'. 

G{n, t\no, 0) = ^p{n, t\no, o). (lo) 


A. Depolarized field autocorrelation function for rotational diffusion 9iXXtiQ,t) 


In a depolarized light scattering experiment, the sample is illuminated by linearly polarized light with the electric 
field vector oriented in z-direction. Both the scattering vector Q and the electric field vector of the scattered light are 
perpendicular to the z- direction. In VH-geometry, scattered light is only detected if the electric field vector changes 
its orientation during the scattering process by a rotation of the scattering objects. To calculate the statistical average 
over the initial orientation Dg and the final orientation D, the transition probability p(r2, t|r2g, 0) between these two 
orientations is needed. For a freely diffusing particle it can be obtained from the rotational part 


£rotp(^^,f|f^0,0) = DroJ^p{^,t\no,0) = 


XL 

m 


p(D,t|r2o,0) 


of the Smoluchowski equation ([^. The formal solution of Eq. (11) is 


( 11 ) 


OO I 

p(D,t|Dg,0) = ^ ^ exp(-D,ottP^Yi*^mYim{^o) 
l—O m— — l 

OO I 

= E E Yiimru^o) ■ ( 12 ) 

l—O rn— — l 
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Combining Eqs. (§, and ([l0|) leads to the expression 


9lUQ,t) = 


X) fhm.iiQ)Yi*mii^o) G{fl,t\no,0) ^ /i 2 m 2 (Q)^/ 2 m 2 (^) 

hmi l2m2 / f2n,f2 


S fhmiiQ)Yi*jni(^o) Pi^o) X] //2m2 (Q)^i2m2 (^o) 

, ilTTli l2m2 / f2o 

§§EfimiQ)Yi*^i^o) E expH(Z + l)i^,ott) E /i"m"(Q)>i"m"(f^)df2dOo 

47r 47r Im I'm' l"m" 


§Efim{Q)Y*^i^o) E //w(Q)i1w(f^o)df^o 

47r Im I'm' 


for the field autocorrelation function. Here, dfJ = sin^?d■dd(^ is the differential for the orientational integrals. 
With the orthogonality relation 


(13) 


Yi*^{n)Yi,^,in)dn = 6u'S„ 


47r 


one obtains 




EE E flfniQ)exp{-l{l + l)Drott)6w5mm'fl"m"iQ)6l'l"6r, 

Im I'm' l"m" 

E E flm{Q)fl'm'{Q)Sll>Smm' 

Im I'm' 


E fLiQ) exp i-l{l + l)£>rott) 

Im 

EfLiQ) 


1 


{PiQ,^))n 


Im 


(14) 


fLiQ) (-1(1+ ^)Dy:ott). (15) 


The orientation-averaged particle form factor 

p{Q) = {PiQ,n))^ = J2fLiQ) 


Im 


(16) 


is the scattered intensity of a randomly aligned ensemble normalized to the forward scattering, i. e., limQ_,,o P{Q) = 1. 


B. Depolarized field autocorrelation function for translational diffusion: the isotropic part 

9l^Dtr) i^^ '') 


By means of the short-time limit of the orientational transition probability density 

limG(D,t|Do,0) = = ^<5(0 - Oo) (17) 

Im 

the contribution of translational diffusion to the field autocorrelation function is accessible. Physically, the limit 
t —)■ 0 means that the orientational degrees of freedom are frozen. The initial condition p(r, 0) = i5(r) for translational 
diffusion reads p(Q, 0) = 1 in Fourier space. Hence, with the notation Q = QD the solution of the partial differential 
equation 


P{D,,)PiQ,^,t) 


-Q'^{Dtr)p{Q,^,t) 


dt 


p{Q,n,t) 


(18) 


can be written as 


p{Q, n, t) = exp (-Q^(Ar)i) , (19) 

since a scattering experiment selectively probes the projection of the diffusive path to the scattering vector Q. 

Using Eq. ( [T7| ), the depolarized field autocorrelation function for the orientation-averaged translational diffusion of 
an ensemble of randomly aligned particles reads 
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exp {-Q'^ {Dtr)t) fl’m'{Q)Yl>m'{^)) 

Im I'm' / < 


Oq ,o 


E fimiQ)Y;jno) p{no) E fi'm'iQ)Yvm'ino) 


( 20 ) 


After performing one orientational integration, the numerator and the denominator in Eq. (20) are identical except 
for the exponential factor. Thus, one obtains 


(Q: t) = exp (-Q^(Ar)i) ■ (21) 

This contribution is the only one also visible in polarized dynamic light scattering experiments, where the polarization 
of both incident and scattered beam is vertical (VV-geometry). As a consequence, such experiments are only sensitive 
to the orientation-averaged translational diffusion coefficient. 


C. Depolarized field autocorrelation function for translational diffusion: the anisotropic part 

The final contribution 5^ADtr(Qi^) depolarized field autocorrelation function originates from the anisotropy 

of the translational diffusion. This contribution has only been addressed for infinitely thin rods so far.^^ Here, we 
provide a general expression in terms of rotational invariants. The partial differential equation 


A(cosi5)p((5, r2,t) = ^p{Q,n,t) 

resulting from the traceless symmetric tensor has the formal solution 

p{Q,fl,t) = exp r-^(5^AArA(cosi?)tj . 


( 22 ) 


(23) 


The second Legendre polynomial can be written as P 2 (cosi?) = (47r/5)^/^l2,o(J^) in terms of spherical harmonics. 

Herewith, the contribution of the depolarized scattering function related to the anisotropy of the translational 
diffusion is given by 


E/im(Q)l)*„(f^o)4^<J(^^-^^o)exp( -^Q2^Ar( ^Vl2.o(W E(Q)Am'(f^) 


i Im 


5 J 


I'm' 


glADjQ^t) = 


Eflm{Q)Y;Jno)p{no) E ^m'(Q)Arn'(^^o) 


■ (24) 


The short-time behavior of g^^Dt iQ’ 0 is obtained by expanding the exponential function in a Taylor series. In the 
short-time limit we obtain 


2 / Att \ 2 

p(Q,D,t) = l--Q^ADtr f yj Y2,omt + Oie). 

With this expansion the depolarized scattering function gi^D^^{Q,t) reads 


(25) 




§EfimiQ)Y;jno) E /iw(Q)A™'(f^o)ciDo 

477 Im I'm' 

§j:fimiQ)Y*^ino) E /iw(0)W(f^o)clDo 

477 Im I'm' 




flm{Q)Y*^{ClQ)Y2fl{flo) E fl'm'{Q)Yl>m'{^o) dfto 

477 Im I'm' 

§Efim(Q)Y;^(^o) E fl'm'(Q)Yl^^'(^o)dno 

477 Im I'm' 


+ Oif ), (26) 
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where the integration with respect to ft has already been performed. Using the orthogonality relation (14), this leads 
to 


dl.ADt, 


iQ,t) = l- 


2Q^ADt,t 

3P{Q) 




Im Vm' 


(Q) <#) Yi*^ino)Y2,oino)fum'{Q)Yi,m'ino) df2o + 


(27) 


477 


The remaining integral over three spherical harmonics depending on fJo can be solved employing Gaunt’s theorem.l^ 
Thus, one finally obtains 


9i.AD,AQ^t) = 1 - 


2Q‘^ADtrt /47r\ 


3P(Q) V 5 y 


EE 


Im I'm' 


' C{1', 2, 1- m', 0, m)C{l\ 2, t, 0, 0,0) + 0{f) 


= 1 - 


2Q^ADt^ t 

3P(Q) 


EE flm{Q)fl'm'{Q) 


Im I'm' 


21 ' + 5 

2Z + 1 


C{1', 2, 1; m, 0, m)C{l', 2, 1; 0,0, 0) + 0{t^) 


(28) 


with the Clebsch-Gordan coefficients C(li,l 2 ,h]‘rni,m 2 ,rn^). 

For solids of revolution with an equatorial mirror plane, only the coefficients /;_o with even I differ from zero. Hence, 
in this case just contributions 


C{Q) = 


2 

3P(Q) 



C^{2l',2,2l;0,0,0)f2i,o{Q)f2i',o{Q) 


(29) 


have to be considered. Only Clebsch-Gordan coefficients (7(2^, 2, 2Z; 0, 0, 0) with \2l' — 2l\ <2 are nonzero. The 
relevant Gaunt coefficients appearing in the coupling function of solids of revolution with an equatorial mirror plane 
are 


ClM 


fil' + l 


1/2 

C^{2l',2,2l-0,0,0). 


With these coefficients the function C(Q) can be written as 


C{Q) 


2 

mQ) 


EE flm iQ) fl'm' (Q) 

Im I'm' 


leading to the expression 

9YadJQ, t) = l- C{Q)Q^AD,,t + 0{e) 


(30) 


(31) 


(32) 


for the short-time behavior of the field autocorrelation function resulting from the anisotropy of the translational 
diffusion. Iteration for a finite time t then yields 

9iADtr (e, t) = exp (-C'(Q)Q^AAr t) . (33) 

Unlike the weighting factors for the rotational contribution [see Eq. (HH] , the coupling function can be calculated via a 
different route avoiding the expansion in spherical harmonic^^ as the orientation average of the form factor weighted 
by 2P2(cos'd)/3. If P{Q,'d) denotes the form factor of a particle with the orientation 11, the coupling function C{Q) 
for solids of revolution can be written as 


CiQ) = I P(Q,i/)2p,(cosi?) sini^dd, (34) 

0 


since the angular dependence of the diffusive quadrupole contribution (Z = 2) just leads to a factor 2P2(cosi?)/3. In 
the rotational part, however, the multipole contributions are weighted by Z(Z -|- 1), which makes a separation of the 
multipole contributions necessary. 
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In the context of scattering from solids of revolution, the term translational-rotational coupling denotes a dependence 
of the translational mobility and the orientation of the particle in the laboratory coordinate system defined by the 
scattering vector: if a prolate particle has the orientation SI at time t, its mobility parallel to the orientation vector 
is larger than perpendicular. Since scattering experiments, however, probe the projection of the diffusive path to the 
scattering vector, also the rotation of anisotropic particles influences the scattered intensity related to translational 
diffusion. 

As a side remark, we note that for objects of more complex shape, where due to the absence of inversion symmetry 
nonzero off-diagonal elem ents i n the diffusion tensor exist, in terms of hydrodynamics another kind of translational- 
rotational coupling occursIn that case, a rotation causes a translation and vice versa, leading to a translational- 
rotational coupling even in the particle’s frame of reference. Such effects, however, are not present for shapes with 
inversion symmetry. 


III. EXPANSION OF THE SCATTERING FUNCTION IN SPHERICAL HARMONICS 


For the calculation of the depolarized scattering functions in Eqs. (151 and (28) the expansion coefficients fim{Q) 
are required. These quantities are Hankel transforms of the expansion coefficients The latter quantities are 

calculated in this section. 

Let us define the reduced scattering length density of a solid of revolution in terms of polar coordinates 




1 : inside 
0 : outside 


(35) 


of the particle. Assuming axial symmetry in z-direction, this function is independent of the azimuthal angle Lp, hence 
= 4'(r, ■!?). Let us additionally assume inversion symmetry related to the particle’s center of mass. Then, 


together with the axial symmetry, an equatorial mirror plane is generated leading to the parity relation 

dt (r, i?) ='I'(r, TT —-d). (36) 

The spherical harmonics Urn(S^) with 

Um(d, 1 /?) = f !! , Pim(cosd)e*’"‘^ Vm>0 and (37) 

\ 47r (/ -|- m)!y 

= Vm<0 (38) 

form a complete orthogonal basis. Hence, the reduced scattering length density can be expanded in spherical harmonics 

oo ra 

5((r,d)=^ X! fim{'r)Yi^{-d,(p) (39) 

m——l 

with the expansion coefficients given by 

fUr) = (jj)Yr^in)^ir,n)dn. (40) 

477 


A. Intersections of solids of revolution with spheres 

With knowledge of the meridian curve x{z) of a solid of revolution, the function 'I'(r, d) can easily be obtained 
(see Fig. [^, if the polar angle dc(r) of the intersection of a sphere with the radius r and the solid of revolution is 
calculated. In dependence on dc(r) the function 'k(r, d) can be written as 

'I'(r, dc) = I — 7t(cosdc) -I- 7t(cos(7r — dc)), (41) 

where 77 (cos d) denotes the Heaviside function 


77(cosdc) 


<5(cos(d — dc)) d cos d. 


-1 


(42) 
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FIG. 1. The reduced scattering length density T(r, i?) as a step function of the critical pole distances and tt —i9c exemplarily 
shown for spherocylinders. The critical pole distances i9c and tt — depending on the radius r are accessible from the meridian 
curve x{z) of the particles. 


For the calculation of the expansion coefficients the pole distances of intersections of spheres and the solids of 
revolution are necessary. Let ^ = cos'dc be the cosine of the critical pole distance where the sphere intersects the 
solid of revolution. Due to the equatorial mirror plane, an identical intersection can be found at tt — 'dc- The critical 
polar angles are accessible via the meridian functions x(z). Let us for simplicity use reduced units, i.e., we set the 
equatorial radius to i?eq = 1- 

Exemplarily, we discuss the expansion for spherocylinders and double cones here, while cylinders, ellipsoids, and 
spindles are addressed in Appendix A. A spherocylinder with the aspect ratio v consists of a cylinder with the radius 
i?cq and the length 2i?eq(j^ — 1) with two hemispheres with centers at {iz — l)i?eq and —{v — l)i?eq (see Fig. [^. 

Hence, the meridian curve for a spherocylinder with the rotation axis in ^-direction is 


x{z) = 




1/2 


: -(J 2 - 1) <z<v-l 
: \z\>v-l 


The symmetrically equivalent pole distances of intersections with spheres can be written as 


(43) 


^ = cos = < 


± 

± 



A (l-(^-l))' 

r2 


1/2 


r< ((i^-1)^ + 1)" 

{{v-lY + lY <r<v 


(44) 


In the expansion of double cones with the meridian curve x{z) = 1 — \z\/v (see Fig.[^, more than two symmetrically 
equivalent intersections with spheres can exist for Tc < r < 1 with = sin(arctan v) denoting the radius of the largest 
sphere completely inside the double cone. In this region, the cosines of the critical pole distances read 


6.2 = COS!?,,,, = ± (^'(l + ■ (45) 

For 1 < r < y only two symmetrically equivalent intersections exist at the critical pole distances 

^ = cosdc = (l + (^"(l + • (46) 


B. Expansion coefficients 


The angular-dependent function 4'(r, (/?) can be expanded by using the expansion coefficients defined in Eq. (401. 


Since for rotationally symmetric solids of revolution, 'I'(r, t?) does not depend on the azimuthal angle (/3, the integration 
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FIG. 2. Meridian curves for spherocylinders (Ihs) and double cones (rhs) as well as critical pole distances for intersections with 
spheres. 


over ip gives a nonzero contribution, namely a factor 27r, only for m = 0. For m ^ Q the azimuthal integral vanishes 
so that only the coefficients fi,Q are relevant. 

The solids of revolution considered here have an equatorial mirror plane, i.e., lE'(r,d) = 'I'(r ,tt — d). Due to the 
parity 

Pi{x) = i-lYm-x) (47) 


of Legendre polynomials only coefficients with even I differ from zero. 
Thus, the relevant expansion coefficients for prolate particles are 


f 2 i{r) = 2TrJ P2i(a;)4'(r,x) da; 


-1 


47r 


-i+? 1 

= ((4Z + l)7r)5 I J P2i{x)dx+ J P2iix) dx 

-1 i-c 

1 

= ((4Z + l)7r)5 2 J P 2 /(a;) dx = 2 ((4Z + l)7r)5 S 2 i(^(r)), 
i-C 


(48) 


where x = cost? and ^ = cosdc. The coefficients S(^(r)) depend on the shape of the solids of revolution only via the 
cosine of the critical polar angles ^(r) = cos dc(r). 

For double cones, in the region Xc < r < 1 additional intersections centered at the equatorial plane occur, leading 
to additional contributions 


^2i{Cii’)) = 27r j P 2 i{x)dx = 2 {{41 +l)Tr)^ J P 2 i{x)dx (49) 

0 

to the expansion coefficients, which read in this case f 2 i{i~) = 2((4^ + l)7r)^/^(E!2i(C(?’)) + A(^'(r))). As a test for the 
expansion in spherical harmonics, the meridional cross sections of the solids of revolution are reconstructed from the 
expansion coefficients fi{r). In Fig. the sum of the contributions up to the order Z = 40 is exemplarily displayed 
for spherocylinders and double cones. 
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FIG. 3. Reconstruction of a spherocylinder (Ihs) and a donble cone (rhs) with an aspect ratio of = 5 from the expansion in 
spherical harmonics up to the order I = 40. In the scattering functions, the deviations from a nniform scattering length density 
are relevant for large scattering vectors probing the short length scales of correlations in the scattering length density. 


IV. SCATTERING FUNCTION 


The static scattering function of a particle with a homogeneous scattering length density can—normalized to the 
forward scattering—in terms of its form factor be written as 


P{Q) = P{Q,n) = 


1 




(50) 


where Q denotes the scattering vector and V the volume of the particle. The function lE'(r) is the reduced scattering 
length density with the values 'l'(r) = 1 inside and lE'(r) = 0 outside of the particle. The normalization to the squared 
particle volume relates the scattering intensity at the scattering vector Q to the forward scattering at Q = 0. 

With 'l'(r) = 'I'(r’r2), where denotes the direction of r, using the expansion in spherical harmonics the form factor 
reads 





I / (^E/2;W>?co(^,<P)e*'^"“^’'dcosdd(^dr 
' ) 

dcosd dipdr 1 . 


(51) 


Here, the direction of Q is chosen to be the positive z-direction. Then the angle enclosed between Q and r is 'd. 
Under this assumption the expansion coefficients f 2 i{Q) in Fourier space can be obtained via Hankel transforms 


f2i{Q) = {if‘ J 47rrV2i(?’)j2i(r)dr , (52) 

0 


where j 2 i{r) denotes the spherical Bessel function of the order 21. Since these expressions do not depend on angular 
coordinates anymore, the average over all orientations is directly obtained, and the result is independent of the choice 
of the laboratory coordinate system. 

The scattered intensity at the scattering vector with the modulus Q and the direction is 

P(Q,II) = ^ (^^f2r(Q)r2i',o(n)^ ■ (53) 
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Hence, for an ensemble of statistically aligned particles with the odf p(r2) = l/(47r), the orientation-averaged form 
factor is given by 




1 


Y2 



4-k 


E E hiiQ)f2AQ)Y2iomY^*',om 

i i' 


E f^liQ) I 

I 


(54) 


since the odf eliminates due to the normalization to the forward scattering of an isotropic ensemble. 

The form factors can be obtained via a different route by integration in cylinder coordinates.^^ The semianalytic 
integration leads to form factors P{Q, i?) depending on the direction of the scattering vector Q and the direction of 
the particle axis f with cosi9 = Q • f. The orientation average P{Q) = {P{Q,'&))■& can be calculated by numerical 
integration. 

Let us define reduced expansion coefficients 


SAQ) 


fhiQ) 

E/KQ) 

I 


fUQ) 

vAP{Q,mn ■ 


(55) 


Then, SoiQ) is formally the reduced monopole contribution to the scattering intensity, 52(Q) the quadrupole contri¬ 
bution, SAQ) the hexadecapole contribution, and so on. 

If we further define 


eiiQ) = J2Sn{Q) 


(56) 


n—0 


as the sum of the first I expansion coefficients, the limit lim;_>oo ©((Q) = 1 is a test for the convergence of the 
expansion if (P(Q,ll))o is calculated via numerical integration in cylinder coordinates. 

In Figs. |^and[^ the first expansion coefficients So{Q), S 2 {Q), and S^Q) are displayed for various particle shapes. 
We focus on cylinders, spherocylinders, spindles, and double cones since the convergence of the expansion of ellipsoids 
with large aspect ratios in terms of spherical harmonics is rather unrewarding. The plots in the lower right corners 
of Figs. 1^ and illustrate the convergence of the expansion exemplarily for spherocylinders by showing different 
Gi{Q). For highly anisotropic objects, the limitation of the Q-range is visible even up to the order I = 30. Expansion 
coefficients of higher order lead to numerical instabilities in the numerical Hankel transform with double precision 
arithmetics. 

The coupling function C{Q) defined in Eq. also depends on the expansion coefficients f 2 iiQ) weighted by 

Gaunt coefficients. Coupling functions for differently shaped objects are displayed in Fig. for the aspect ratios 
V = 2 and v = 6. The increasing similarity of C{Q) for cylinders and spherocylinders with increasing aspect ratio 
is reasonable since the differences in the shape vanish in this case. Only at large scattering vectors probing short 
intraparticle correlation lengths, significant differences exist. 

Combining Eqs. (15), ( [^ , ( [^ , and 0 leads to the expression 


gi^iQ,t) = exp (-Q^ (Ar)t) exp {-C{Q)Q‘^ADtrt)^Si{Q)exp{-l{l + l)D^ott) 


(57) 


for the total depolarized field autocorrelation function since, due to the rotational symmetry of the particles only 
expansion coefficients Si{Q) = fA{Q)/{V‘^{P{Q,^)}n) are nonzero. 

Experimental access to the short-time behavior of the intermediate scattering function gY^{Q,t) is enabled by its 
first cumulants r((5), i.e., the negative initial slope with respect to t. 


r(Q) = - liin 

t—J-O 


dgYAQA 

dt 


= — lim 

t —^0 


d\ngY^{Q,t) 


— Q‘^{P)tY + ADtrC{Q) + 


dt 
III -f 


l)D,otSi{Q). 


(58) 


While in a polarized scattering experiment in VV-geometry the first cumulants are proportional to according to 
the Landau-Placzek relation,in a depolarized scattering geometry two additional contributions exist. On the one 
hand, due to the anisotropy of the translational diffusion, an additional contribution proportional to which is 
however modulated with the coupling function C{Q), occurs. On the other hand, the first cumulants are influenced 
by the rotational diffusion. 
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FIG. 4. Expansion coefficients S 2 i(Q) for prolate solids of revolution with an aspect ratio of = 2.0 and a convergence test for 
spherocylinders with the aspect ratio v = 2.0 (lower right corner). QR^q is the reduced scattering vector with f?eq denoting 
the equatorial radius of the particles. 


V. CALCULATION OF THE DIFFUSION COEFFICIENTS 


The intermediate scattering function f) [see Eq. (57)] as well as its first cumulants r((5) [see Eq. (58)] depend 

on the translational and rotational diffusion coefficients of the considered particles. We calculate these coefficients 
for various particle shapes based on a bead model using the software This software has originally 

been designed for the prediction of hydrodynamic properties of macromoleculesf^but it can also be used for colloidal 
particles^ modeled by a large number of beads.l^ To calculate the diffusion coefficients of such particles, first the 
bead positions have to be determined. In order to achieve a high accuracy, it is important that the bead model 
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FIG. 5. Expansion coefficients S 2 i(Q) for prolate solids of revolution with an aspect ratio of = 6.0 and a convergence test for 
spherocylinders with the aspect ratio v = 6.0 (lower right corner). 


closely represents the specific particle shape. As the number of beads that can be used for the calculations with 
the software HYDRO++ is limited, it is most efficient to place the beads only on the surface of the p article to be 
modeled. This method is referred to as shell model and was originally introduced by Bloomfield.lS2E2l enables a 
much better resolution of the particle shapes than filling models with a full-volume bead representation as beads with 
a significantly smaller radius can be used. Moreover, in many situations shell models even provide more accurate 
results than filling models with the same bead size.^^ Best results are obtained if the beads are placed on the surface 
without larger voids so that they ideally touch each other. However, it is important that the beads do not overlap 
because the software may provide wrong results in that case. If bead representations with overl apping spheres are 
required, the program HYDROSUB can be used to perform the hydrodynamic calculations! ^^ l ^^ l This software also 
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FIG. 6. Coupling functions C(Q) for cylinders, spherocylinders, spindles, and double cones with aspect ratios i/ = 2.0 (Ihs) 
and ty = 6.0 (rhs). 



FIG. 7. Bead model representations of the considered prolate particle shapes with aspect ratio v = 2: (a) cylinder, (b) 
spherocylinder, (c) ellipsoid, (d) spindle, and (e) double cone. 


allows for the implementation of ellipsoidal and cylindrical subunits.l^ In order to include the effect of a wall or a 
substrate, it is in principle possible to modify the calculations by using the Stokeslet close to a no-slip boundary.l^ 
This is in particular relevant to experimental situations with a quasi-two-dimensional setup.l^^^^ 

As we are dealing with rotationally symmetric particles, we apply the following efficient method to calculate the 
bead positions that are required as an input for HYDRO-1— 1-.^ The beads are placed on several rings centered on the 
z-axis, which coincides with the symmetry axis of the particles. These rings have a distance of 2r each, where r is 
the radius of a single bead. We always use the same value of r for all beads representing a particle. The radius R{z) 
of the rings determining the bead positions is given by the meridian curves presented in Sec. Ill and in Appendix A, 
respectively. For cylinders, it is necessary to add beads on the top and bottom areas, as otherwise beads are only 
placed on the lateral surface area of the cylinder. Therefore, additional rings with radii equidistantly distributed 
between zero and the equatorial radius i?eq are considered in this case. 

Once the radii of the rings are known, the explicit coordinates of the beads can be calculated. The maximum 
number of beads that can be placed on a ring with radius R is determined by means of the formula s = 2Rsin{a/2) 
for the chord length s. Here, a is the corresponding central angle. Using the condition s = 2r one obtains that the 
minimal ring radius 


i?min(Af) = W 


sin(7r/A^) 


(59) 
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FIG. 8. Translational (Ihs) and rotational (rhs) diffusion coefficients for cylinders, spherocylinders, ellipsoids, spindles, and 
double cones with i?oq = 50 nm in dependence on the aspect ratio v. The dashed lines represent the analytically derived 
diffusion coefficients for ellipsoids.^ 


is required in order to place iV beads of radius r on the ring. Based on Eq. (59), for each value of R(z) the maximum 
number of beads can easily be determined. Using polar coordinates, the obtained number of beads are positioned 
equidistantly on the ring. The resulting bead model representations for the various considered particle shapes are 
visualized in Fig. ,while the obtained values for the translational and rotational diffusion coefScients are displayed 
in Fig. [8| 


VI. DISCUSSION OF THE FIRST CUMULANTS 


Using the diffusion coefficients shown in Fig. the cumulants for polarized and depolarized scattering experiments 
can be quantitatively analyzed. They are displayed in Fig. for particles with an aspect ratio oi v = 2 and i/ = 6, re¬ 
spectively. For objects with typical lengths comparable or larger than the wavelength, also the rotational contribution 
depends on the scattering vector due to the formal multipole expansion of the scattering intensity. 

The first cumulants as the initial slope of the field autocorrelation function provide the dynamic quantity which can 
be accessed with the highest accuracy from experimental data within the regime of the lowest signal to noise ratio. 
Since the coupling function is calculated as a short time expansion, in addition, the initial slope of the autocorrelation 
function is well described by the theoretical approach discussed in this paper. 

where the first cumulants r^''((5) are compared to the 


The influence of the coupling function can be seen in Fig. 


10 


results neglecting the coupling function. The latter case is represented by the dashed lines. For spindles and double 
cones as objects with pronounced apices, the influence of the translational-rotational coupling cannot be neglected at 
large scattering vectors. If the coupling function C{Q) is negative, for prolate objects with AUtr > 0 the quadrupole 
contribution of the diffusion tensor leads to a slower relaxation of the depolarized scattering function. For cylinders 
and spherocylinders strong modulations of C{Q) occur which lead to an accelerated relaxation if C{Q) > 0. 


VII. CONCLUSIONS 

The macroscopic rheological properties of fluids are considerably influenced by the mesoscale dynamics of suspended 
particles. For shape-anisotropic particles, in addition to translational diffusion, rotational diffusion is relevant for the 
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FIG. 9. First cumulants for differently shaped particles with aspect ratios of = 2 (Ihs) and v = Q (rhs). The first cumulants 
of depolarized scattering r'^^((5) are compared to those resulting from polarized scattering r^'^((3). 


rheological behavior, which depends not only on the volume fraction but also on the aspect ratio and shape of the 
particles. 

Depolarized dynamic light scattering gives access to both translational and rotational diffusion of suspended par¬ 
ticles. Even in highly dilute suspensions without particle interactions, the depolarized scattering function shows a 
nontrivial dependence on the scattering vector Q if at least one typical particle dimension is comparable to or larger 
than the wavelength. The first cumulants of the intermediate scattering function deviate from the Landau-Placzek 
behavior with a proportionality to . The modulation of the cumulants strongly depends on the particle shape and 
is related both to the rotational diffusion and the translation-rotation coupling. The contribution of the coupling 
to the initial relaxation rates can either be positive or negative, depending on the particle geometry and the wave 
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FIG. 10. First cumulants for differently shaped particles with aspect ratios of i/ = 2 (Ihs) and v — G (rhs). The influence of the 
coupling function is shown by comparison of the first cumulants r'^^((5) resulting from depolarized scattering with coupling 
function C{Q) and neglecting the coupling function (C(Q) = 0). 


vector. For prolate objects, AH = Hy — D± is positive. For the calculation of gY^{Q^t) a multipole expansion of 
the scattering function is required. With increasing aspect ratio of the particles an increasing number of expansion 
coefficients is needed to achieve a satisfactory convergence for an appropriate range of scattering vectors. However, 
in light scattering experiments, the maximum accessible scattering vector is limited to Qmax = 47rn/A, where A is 
the wavelength and n the refractive index of the suspending medium. With an effective wavelength of A/n « 400 nm 
for a green laser in aqueous suspensions of par ticles with R^q = 50 nm, which is approximately the typical equatorial 
radius of spindle-shaped hematite particles ,^23 the maximum reduced scattering vector is roughly QRcq « 1.5. Both 
translational and rotational diffusion coefficients of spindles and double cones are due to the smaller surface and 
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volume of these objects significantly larger than those of spherocylinders and cylinders. As a consequence, also the 
difference AZ^tr is larger for the former shapes with pronounced apices. Hence, the effect of the translation-rotation 
coupling visible in the first cumulants is for these particles larger than for the other ones. Hematite particles are 
a promising model system to investigate the coupling between rotational and translational diffusion. However, the 
size of these particles is too small to cover a sufficient Q-range employing visible light as a probe. With coherent 
and linearly polarized X-rays from synchrotron sources and upcoming free electron laser sources, however, practically 
no limitation for the maximum scattering vector in the Q-range relevant for mesoscale objects exists. Employing 
X-ray photon correlation spectroscopy with polarization analysis,^ a sufficiently large range of scattering vectors can 
be accessed to address both translational and rotational diffusion of anisotropic particles. From the first cumulants 
determined over a sufficiently large range of scattering vectors, the quantities (ZZtr), AZZtr, and D-^ot relevant for the 
diffusion of anisotropic particles can be obtained experimentally. 

Our approach can easily be extended to core-shell particles as demonstrated in Appendix B. In the case where 
an anisotropic core is embedded in a spherical shell, also systems with isotropic topology and anisotropic scattering 
power can be treated. Since for a spherical shell AZZ = — D± vanishes, a rotation-translation coupling is not 

present for such a system. Only the multipole expansion of the anisotropic scattering power modulates in this case 
the Q^-proportionality of the first cumulants. 

The approach can also be adapted to describe shape-anisotropic particles consisting of optically anisotropic materials 
with a tensorial polarizability. As long as the diffusion tensor is diagonal, also scattering functions for shapes without 
rotation and inversion symmetry can be calculated. In this case, additional expansion coefficients have to be 

considered. Without optical inversion symmetry, also coefficients with odd order of I are nonzero. 

Moreover, for long rods, effects of particle flexibility may get important. Fluctuations in the particle shape can be 
incorporated into the present analysis by considering a smeared averaged static form factor. The full description of 
flexibility effects beyond this simple effective picture, however, requires a nontrivial extension of our theory which is 
left for future studies. 

Our analysis provides a systematic framework to interpret depolarized light scattering data with respect to details 
in the shape of prolate particles. This is of great importance to the precise characterization of suspensions (beyond 
their aspect ratio) which contain anisotropic nano- or colloidal particles. One important example are ferrofluids which 
involve anisotropic magnetic particles in a carrier fluid and are of great importance to many magnetorheological 
applications.!^ 

Future work needs to be performed for oblate particles and for particles of arbitrary shape that are not rotationally 
symmetric. The latter case requires a more sophisticated desc riptio n for the orientations (e.g., by using Eulerian 
angles) such that the Brownian motion is much more complicated.^SSESl Another line of future research is to incorporate 
the correlations at finite particle density. Particularly interesting is the high-density regime where the suspension is 
getting crystalling^ or glassy. An explicit computation of the scattering function will, however, become much more 
complicated than the analysis presented here due to the intricate character of the orientation-dependent interactions. 
Still, in phases with strict periodic order, such as fully crystalline phases or plastic crystals, one can hope that the 
analysis is somewhat easier since one can apply natural periodic boundary conditions. An important question is 
whether one can achieve a quantitative level of understanding of the role of translational and rotational fluctuations 
of the particles in the stability of specific ordered phases. Finally, an external field can be applied which couples to the 
orientational degrees of freedom and gives rise to new phase phenomena, such as a paranematic-nematic transition,!^ 
new kinds of crystals,!^ and an introduction or destruction of biaxial orientational order (see, e.g.. Refs. 1751 and [75ll . 
An extension of our analysis to the case of noninteracting particles in an orienting external field should be feasible 
and will be a topic of our future research. 


Appendix A: Meridian curves for solids of revolution 


The 
with a 


meridian curve for a cylinder simply is a:(z) = 1. For r < v only one symmetrically independent intersection 
sphere exists, for v < r < 1 an additional circle exists as intersection. The critical pole distances are 



r < V 


V < r < y/v"^ -\-1. 


(Al) 


19 













FIG. 11. Reconstruction of the meridional cross section of a core-shell structure consisting of a spindle-shaped core with a 
length L = 0.75(Ts and an aspect ratio of = 3 in a spherical shell with the outer diameter Gs- The scattering length density 
of the core is pc = 1 and that of the shell ps = 0.2. 


In the region v < r < + 1 the two symmetrically independent intersections lead to the expansion coefhcients 


f 2 i{r) = 2y/{Al + 1) IT - E 2 i{^ 2 {r))]. 


(A2) 


Here, |Ci(^)l < l? 2 (?’)| denote the cosines of the two independent critical pole distances. 
The meridian curve of an ellipsoid is given by 



Only one symmetrically independent intersection with a sphere at the critical pole distance 

- 1 \ 5 

e = eosr?e=(^^j 


(A3) 


(A4) 


exists. 

A spindle with the equatorial radius i?eq and the length 2vRcq can be described as a solid of revolution resulting 
from the intersection of two circles with the radius Rcq{v‘^ + l)/2, the centers of which have a distance Rcq{v^ — 1). 
For the volume the analytical expression 



is obtained.^ The corresponding meridian curve 


= 1 - + I + 2) - 4z2 + 1) 5 

leads to one symmetrically independent critical pole distance 


^ = cos = ± 


/ 

\r (1 — 1 / 2 ) 



(A6) 


(A7) 
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FIG. 12. Expansion coefficients S 2 i{Q) for core-shell particles consisting of a spindle with the aspect ratio i/ = 3 embedded 
in a spherical shell for different lengths L/as of the spindles. In the lower right corner, exemplarily a convergence test of the 
expansion for L/as = 1 is displayed. 


Appendix B: Core-shell particles 

The approach of expanding form factors in terms of spherical harmonics can be extended to core-shell structures. In 
the special case where an anisotropic core is embedded into a spherical shell, optically anisotropic structures undergoing 
isotropic diffusion processes can be described. The translational and rotational diffusion coefficients determined by 
the spherical shape of the shell can be calculated employing the Stokes-Einstein and Stokes-Einstein-Debye relations. 
Since the translational diffusion tensor is spherically symmetric, with AD = 0 the coupling function is irrelevant for 
such a core-shell structure. The formal multipole expansion of the anisotropic core’s scattering power, however, leads 
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FIG. 13. First cumulants for core-shell particles consisting of a spindle-shaped core with an aspect ratio oi i/ = 3 and a 
scattering length density of = 1 and a spherical shell with an hydrodynamic radius of Rs = 200 nm and a scattering length 
density of Qs = 0.2 for different sizes of the anisotropic core. The relative size of the core is quantihed by the ratio of the length 
L of the spindle and the hydrodynamic diameter as of the spherical shell. 


to a modulation of the Q^-dependence of the first cumulants. 

In Fig. im the meridional cross section of a core-shell structure reconstructed from the expansion in spherical 
harmonics is displayed. Here a spindle-shaped core with the scattering length density pc = 1 and an aspect ratio 
of zz = 3 is embedded in a spherical shell with a scattering length density ps = 0.2. The length of the spindle is 
L = 2z/i?eq = O.SiTs, where ag = 2Rs denotes the outer diameter of the spherical shell and R^q the equatorial radius 


of the core. 


The expansion coefficients for core-shell systems with differently sized anisotropic cores are displayed in Fig. 12 


The smaller the size of the anisotropic core is compared to the spherical shell, the less important are the multipole 
contributions to the scattering power. Although the coupling function is irrelevant for these systems, depending on 
the relative size of the anisotropic core, in the first cumulants deviations from the proportionality to occur, which 


originate from the (^“dependent multipole contributions to the scattering power. In Fig. 13 the first cumulants 


for core-shell particles with different sizes of the anisotropic core are displayed. The radius of the spherical shell is 
assumed to be Rg = 200 nm. The translational and rotational diffusion coefficients are calculated by means of the 
Stokes-Einstein and Stokes-Einstein-Debye relations using the viscosity rj = 8.9 x 10“^ Pas of water at 298K. 
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